Pseudomagnetic fields and ballistic transport in a suspended graphene sheet 
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We study a suspended graphene sheet subject to the electric field of a gate underneath. We 
compute the elastic deformation of the sheet and the corresponding effective gauge field, which 
modifies the electronic transport. In a clean system the two-terminal conductance of the sample is 
reduced below the ballistic limit and is almost totally suppressed at low carrier concentrations in 
samples under tension. Residual disorder restores a small finite conductivity. 
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Introduction. Graphene layers which are one or a few 
carbon atoms thick [ij, \^ combine novel electronic en- 



ergy spectrum ( "massless Dirac fcrmions' 
usual structural and mechanical properties 




and un- 

MM 

(for a general review, see Refs. [10|, [ll|, |12[ ). Originally, 
graphene sheets lying on quartz substrate have been pre- 
pared and investigated but later it turned out that the 
freely hanged membranes of macroscopically large sizes 
can be derived (5|. The electronic transport and me- 
chanical properties of these membranes are being inten- 
sively studied [i, 0, i, i, E, [11, [H, m. They demon- 
strate, in particular, a much higher electron mobility, at 
least, at low temperatures, than graphene sheets on a 



substrate 13|, [iJ, ll6| and extraordinary mechanical stiff- 
ness [8|, which makes them especially interesting for ap- 
plications. However, peculiarities of the freely hanged 
graphene membranes are still poorly understood theo- 
retically. 

In this paper we demonstrate that unavoidable defor- 
mations of the membranes by an applied electric field can 
strongly affect their transport properties. The system we 
consider is sketched in Fig.[lja). The calculated deforma- 
tion ho and conductance G are shown in Figs. [2] and [H 
The left panel of Fig. [2] shows the maximum deforma- 
tion as function of carrier density n and slack AL in the 




FIG. 1: (Color online), (a) Sketch of the model of a sus- 
pended graphene sheet under consideration, (b) Fermi circles 
positions in the Brillouin zone in the leads (left) and in the 
suspended region (right). 





3x10" 6x10" 



FIG. 2: (Color online) Left: Maximum height, ho, for a sus- 
pended sheet of length I = 1/xm as function of carrier density 
for different slacks: AL = 2nm; pentagons (blue), AL = 0, 
triangles (red), and AL = — 2nm, diamonds (green). Right: 
Effective gauge field for the same three values of slack. The 
full magenta curve gives the Fermi wavevector, kp. 



sheet. We define slack as the difference of the equilibrium 
length of the suspended sheet and the distance between 
the clamped ends. Negative AL implies that the sheet 
is under tension (see below). The right panel of Fig. [2] 
gives G as function of carrier density for the same three 
values of slack. In the presence of a finite deformation 
parametrized by the maximum vertical displacement ho, 
the conductance is reduced. We found that at not too 
small carrier concentrations n, function G(n) is given by 
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where W is the width of the sample, kp = ■\/7r|n| is the 
Fermi wavevector, and Ay of dimension of inverse length 
is related to the deformation [Eq. ([4])]. For graphene 
initially under tension [9|, where /iq is nearly constant, 
and low carrier concentrations, we find a nearly complete 
suppression of transport. The physical mechanism of this 
phenomenon can be understood as follows. The deforma- 
tion shifts the Dirac points by the amount Ay [Fig.[IJb)], 
which creates a mismatch between the graphene leads 
and the suspended region. If this shift exceeds the di- 
ameter 2fci? of the Fermi circle, the electrons are fully 
reflected. 

Recent experiments J16l | show a qualitative agreement 
with our Fig. [3l however, at very low n the conductivity 




FIG. 3: (Color online) Conductance as function of carrier 
density. The values of AL and the corresponding symbols are 
the same as in Fig. 2. The full magenta line gives the ballistic 
conductance in the absence of deformation. The width of the 
sheet is W = 1 nm. 



saturates at at a /inite value ^ e^//i. We attribute this 
to residual disorder-assisted tunneling processes. 

The model. We consider a strip of graphene clamped at 
two parallel edges x = ±L/2, see Fig.[T] We assume that 
the width of the strip in the undeformed state is L + AL, 
where AL can be of either sign. The strip is suspended 
above a control gate, whose electric field induces electron 
concentration n in graphene and exerts on it the pressure 
P = {2Tre^n^)/e. 

The resultant height profile h{x) of the strip has a 
well known analytic form 17[ under the condition Hq = 
max/i(x) <S^ L where the linear elasticity theory ap- 
plies. We focus on the case where n is either compa- 
rable or much larger than uq ~ ^(16/7r)(e«;)/(e-^L3) « 
(6 X 109cm-2)/(L/l/xm)3/2, where k k 1.1 eV is the 
bending rigidity of graphene [l8|. In this case the de- 
formation is nearly parabolic: 
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where u = {n/nQj^L/ho)^''^ ^ 1. The maximum defor- 
mation Hq is the positive root of the cubic equation 
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where E w 22eV/A is the Young's modulus of 
graphene [8[. The values of ho for different values of 
the slack are plotted in Fig. [2l 

The deformation of the strip induces perturbations of 
two types acting on electrons: the scalar potential V{x) 
and the effective vector potential [19|, |20| A{x). We ex- 
amine these potentials below, starting with A(a::). 

Vector potential. The aforementioned shift of the Dirac 
points [Fig. [ijb)] is equivalent to the effect of a constant 
vector potential A = {Ax,Ay). We use this latter formal- 
ism in the following as it can be easily generalized to more 
complex situations. The role of A is the largest when the 
"zigzag" direction is along the y-axis. Assuming this is 



the case, we obtain 



A,ix)=0, Ayix)^CiC 
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where /3 = dlog(7o)/dlog(a) « 2 is the dimensionless 
electron-phonon coupling parameter, 70 « 3 eV is the 
nearest neighbor hopping, a « 1.4 A is the distance be- 
tween nearest carbon atoms, ^ = ±1 is the valley index, 
t ~ PL^/(8ft-o) is the horizontal component of tension 
per unit length at the edges, and Ci is a parameter of 
order unity which depends on the relative displacements 
within the unit cell, determined by the microscopic force 
constants [21[. (Note that u = L^Jt/in.) The corre- 
sponding effective magnetic field B{x) = —dxA^ consists 
of two narrow spikes at the edges x = ±L/2 J22J . 

Transport. To compute the two-terminal conductance 
G through a graphene sheet we assume perfect semi- 
infinite graphene leads of the same chemical potential at 
both ends of the strip. Since the perturbations depend 
only on x, the ky momentum is conserved. However, the 
effective magnetic field at the edges shifts the mechanical 
momentum, ky ^ ky — AySig'n{x) , of electrons that enter 
the strip. This leads to a nearly complete reflection of 
electrons at low density where kp < Ay/2. Similar effect 
has been previously examined in the context of transmis- 
sion of Dirac particles through a region of homogeneous 
magnetic field |23j . However, the present problem has 
new qualitative features, see below. 

For a constant Ay^ the transmission coefficient T{ky) 
can be computed analytically: 
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{ky-q)^. If k{Ay)^ < 0, 



then k{Ay) is pure imaginary, so that sin [k{Ay)L] — s- 
— sinh \k{Ay)L\. In this case T{ky) is exponentially 
small. The plot of T(ky) is shown in Fig. H] using the 
= fci^-sin^ for -7r/2 < 6 < 7r/2 and 
The transmission is indeed almost 



parametrization h 



10 



11 , 



zero for a range of incident angles 9. In addition, we see 
Fabry-Pcrot resonances because of the multiple scatter- 
ing off the two interfaces. (Such resonances are essen- 
tially absent if the field is uniform [23[.) 

Neglecting the contribution of edge channels, which is 
permissible when the number of bulk channels kpW/ir is 
large, the conductance can be computed from 



4e2 f dky , ^ 
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This integral can be done analytically in the limit kp 3> 
\Ay\ and Lkp > 1, see Eq. [H Using Eqs. (3)-(6), we 
also computed G(n) numerically for three AL shown in 
Fig. 3. For AL > 0, where the tension and therefore 
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FIG. 4: (Color online). Angular dependence of the transmis- 
sion for ho = 20 nm, no slack, and L = 1 /xm. The carrier 
density is n = 2 x 10^^ cm^^. The two curves correspond to 
the two inequivalent Dirac points. 



I^j^l rapidly increase with n (Fig. 1), the deviations from 
the ballistic formula occur at large n. In contrast, for 
AL = — 2nm, where the tension is approximately con- 
stant, the largest effect of the gauge field is felt at low 
concentrations, leading to a nearly complete vanishing of 
G at \n\ < 2.0 x lO^Ocm-^. 

Disorder effects. Any realistic system contains dis- 
order. When weak, it does not change qualitatively 
the above results for G{n) at high carrier concentration 
\n\ :$> Ay. In the opposite limit, when G is strongly 
suppressed, the effect of disorder is important because it 
relaxes the constraint of momentum conservation, com- 
pensating for the momentum shift due to the gauge field 
at the interfaces x = ±L/2. 

Consider the experimentally relevant case [16| where 
the elastic mean-free path I is comparable to the system 
size L. Then, at 1 < G/{Ae'^lh) < kpW we have a 
regime where the disorder is weak enough that its effect 
on the average conductance is still negligible yet it is 
strong enough to fully mix the transverse modes in the 
suspended region. In this case the conductance is limited 
by the two interfaces, which act as classical resistors in 
series [2J], G = Gi/2. Let us compute the conductance 
of a single interface G,; with and without including the 
effect of disorder. 

If disorder near the interface is neglected, Gi is given 
by the expression similar to Eq. ([6|) where T{ky) is now 
the transmission coefficient through a single interface: 
T{ky) = [4fc(0)/c(A)]/[(fc(0) + fc(A))2+A^]. This formula 
holds if fcy-mode is propagating, 3mfc(0) — 3mfc(Ay) = 
0. Otherwise, it is evanescent and T{ky) = 0. Each 
evanescent mode decays exponentially either to the left 
or to the right of x = interface, depending on which 
of Smfc(O) and 3mfc(Aj,) is nonzero. At kp < \Ay\/2 all 
modes at the Fermi energy are evanescent. Therefore, if 
the disorder is neglected, Gi vanishes. 

Let us now include disorder-induced mixing among the 
evanescent and propagating modes, which gives a cor- 
rection AT{ky) to each T(ky). To the lowest order in 
the concentration Us of scatterers, this correction can be 



Ay+kp oo 



AT(fc,y)= / -^ j dxns5T{ky,k'y,x), (7) 

Ay — kp — oo 



where ST{ky,k'x) is the off-diagonal transmission co- 
efficient due to a single scatterer at position r = 
{x,y). Function AT(ky) has the dimension of length, 
similar to the transport cross-section Ag = l/{nsl). 
On physical grounds, we expect 6T(ky,k'y,x) ~ 
(As/fci?) exp (— 2|fc(Ay)a;|). Here the exponential repre- 
sents the probability of the evanescent wave to reach the 
scatterer and the prefactor provides the correct units and 
scaling with disorder strength. The dominant contribu- 
tion to AT{ky) comes from the scatterers located in the 
strip \x\ < l/\k{Ay)\. Integrating over x, ky, and k' we 
finally get 
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where Ci is a numerical coefficient and o'(n) is the 
conductivity of a sample with size L ^ I. A formal 
derivation based on the Green's function formalism yields 
C2 = 4. For estimate, we can take n = lO^cm"^, 
a - 4eV^, 1^ = l^m, and Ay = 2 x 10^ cm'^ We 
then find G ^ 2e^//i, i.e., an appreciably large value. 

Scalar potential. The deformation of the graphene 
strip also creates a scalar potential V{x) in the system. 
Our estimates below indicate that it is relatively small, 
so that its presence is not expected to change the results 
shown in Figs. [3] and [5] in a major way. We will discuss 
it briefly, for completeness. 

The bare potential induced by the deformation is: 
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where the first term is the change in electrostatic poten- 
tial due to the change in distance to the gate, the second 
term gives the deformation potential induced by a local 



compression |2ll . |25| u^x 



and Vb« 10 eV. 



Within the linear screening theory, the Fourier trans- 
form y of y is given by V{q) = Vextiq) / ^{q) i where e[q) 
is the dielectric function. For reasonable carrier concen- 
trations, the potential in Eq. ^ is smooth over distances 
~ kp"^ -C L. We can use the Thomas-Fermi (TF) ap- 
proximation, e{q) = 1 -|- kg/ \q\, where k^ = Aakp is the 
inverse TF screening length and a = e'^/ehv ~ 1 is the 
dimensionless strength of the Coulomb interaction. The 
screened potential can be computed analytically in terms 
of special functions. In the limit ksL ^ 1 and at dis- 



tanccs greater than k^ ^ from the boundaries it reads: 
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The potential at the edge is given by 
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where C3 ~ 1. Thus, the divergences in Eq. p^ are cut 
off at the distance of the order of the screening length 
1/ks from the edge, as expected. 

For n ~ 3 X 10^^ cm~2, the deformation and the elec- 
trostatic potentials at the edges are comparable in mag- 
nitude and together amount to about 10% of the Fermi 
energy. Although this is not a negligible amount, it is still 
numerically small, so that the linear-response screening 
approach is justified. 

Experimental implications and future directions. In 
this paper we focused on a simplest geometry of a sus- 
pended graphene strip (Fig. [T|) for which the calculation 
of the transport properties can be done semi-analytically. 
We have found that reasonable values of the tension, 
AL/L ~ 0.2% in a suspended graphene strip can lead 
to the almost total reflection for a significant range of 
incoming momenta, which causes a downward shift of 
the two-terminal conductance compared to the ballistic 
limit, Eq. ([1]). Residual disorder can partially compen- 
sate for this shift, which may be the case in current ex- 
periments [13|, [ij, |l6|. Note however, that our theory 
cannot be directly compared with experiments, as the 
sample geometry and therefore the configuration of the 
gauge field can be more complicated than what we have 
assumed here. 

There are a number of possible directions for future 
study. One interesting problem is how the deformation 
would affect the quantum Hall effect (QHE) in the sus- 
pended graphene. Below we offer a preliminary discus- 
sion of this question. 

For the model considered, the effects of the deforma- 
tion are restricted to the x = ±L/2 interfaces. The Lan- 
dau levels near such lines will be modified (exce pt f or the 
A^ = one, which is topologically protected [26|, \2l\). 
Quasiclassicaly, the reflection at the interfaces creates 
skipping orbits, which propagate parallel to the y-axis 
but in opposite directions on the two sides of each in- 
terface. This could lead to backscattering of the edge 
currents and modification of the QHE. The effect is the 
strongest when kp < \Ay\. For low- lying Landau levels, 
where kp ~ 1/^s, an estimate of the external magnetic 
field i?« below which the QHE is affected can be derived 
from the condition Ib{B^,) ~ |Aj,|~^. For AL = — 2nm, 
and Ay^2y. 10^ cm-\ it yields B^ - 0.7 T. 



In a more realistic geometry, a small fictitious mag- 
netic field will also exist inside the suspended region. Its 
magnitude is of the order of 0.05 T for the same AL and 
L. Landau levels mixing in the bulk occurs when the 
external field is comparable or smaller than this value. 
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